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Time-dependent Dalitz plot analysis of B° — ^ D^K^n^ decays. 
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We present for the first time a measurement of tiie weak phase 2/3+7 obtained from a time-dependent 
Dalitz plot analysis of B° D^K°n^ decays. Using a sample of approximately 347 x 10® BE 
pairs collected by the BABAR detector at the PEP-II asymmetric-energy storage rings, we obtain 
2/3 + 7 = (83 ± 53 ± 20)° and (263 ± 53 ± 20)° assuming the ratio r of the 6 -» w and 6 -» c decay 
amplitudes to be 0.3. The magnitudes and phases for the resonances associated with the 6 — »■ c 
transitions are also extracted from the fit. 
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PACS numbers: 



The weak phase 7 = arg ; 



where Vij are ele- 
ments of the Cabibbo-Kobayashi-Maskawa quark-mixing 
matrix [l|, is the least constrained angle of the unitarity 
triangle [2|- Over the past few years, several methods [3[ 
have been employed to measure 7 directly in charged 
B decays 4;], where sensitivity to the weak 

phase arises from interference between the b ^ c (fa- 
vored) and b u (suppressed) transitions. In addi- 
tion, decays to two-body final states containing charm 
have been studied, such as D^*'>^Tr^ and B^ — > 

D^p^ 5] which are sensitive to the weak phase 2f3 + j 

due to B'^-B^ mixing. The phase /? = arg ^— y^y.^ ^ is 
well measured in neutral B decays to charmonium final 
states 0. The sensitivity of this method is limited by 
the ratio r between the b u and b ~f c transitions, 
which is expected to be very small (^ 0.02). Three-body 
B decays have been suggested [tJ as a way to avoid this 
limitation, since r in these decays could be as large as 0.4 
in some regions of the Dalitz plot. 

In this paper we report on the first measurement of the 
weak phase 2/3-1-7 obtained from a time-dependent Dalitz 
plot analysis of the decay D^K'^tt'^ Q (charge 

conjugation is implied throughout the paper). In the de- 
cay i?° D^K^TT^ the three body final state is reached 

predominantly through intermediate B'^—^D Kg and 

B^^D^ K*^ decays. In the first case, D indicates a 
£»S(2400) or a D{{2mQ) state produced through b ^ u 
and b c color-suppressed transitions. In the second 
case, i^*(892), K^{U30), i^2*(1430) and i^*(1680) are 
produced through b ^ c tree-level transitions. A small 
contribution from the b ~* u decay i3''^-D*+(2573)7r~ is 
also expected. 

Defining x as the vector of the two invariant masses 
squared w?{K^tt^) and m^{D^TT^), the amplitude A at 
each point x of the Dalitz plot can be parameterized as 
a coherent sum of two-body decay matrix elements ac- 
cording to the isobar model Q : 



E 



aje'^'BWix;Mj,Tj,Sj), (1) 



where c (u) indicates the 6 — > c (& — s- u) transition and 
(j) is the total strong phase. Each resonance j is param- 
eterized by a magnitude phase Sj, and a factor 
BW{x; Mj,Tj, Sj) giving the Lorentz invariant expres- 
sion for the matrix element of the resonance as a function 
of the position x, the spin s, the mass M, and the decay 
width r. 

The time-dependent probability of a or initial 
state to decay to a final state with a _D+ or can be 



expressed as ; 

Pix,At,trj) = X X 

{1 - rj^Cix) cos(AmdAt) -I- ^S^{x) sin(AmdAt)}. (2) 



Here: 



S,j{x) 
Cix) 



2Im(A,(x)A^(x)e'(^^+T)+'"('^'=(^)-'^"(^») 

A,(f)2-H^„(f)2' 



(3) 



where At is the difference in proper decay times of the 
reconstructed meson Bj-ec and the flavor-tagging meson 
Stag, ^ = +!(-!) if the flavor of the B^ec is a B°{B°) 
and 77 = +1(— 1) if the final state contains a D^{D^). 
We use the world averages for the lifetime tb and the 
mass-eigenstate difference Am^ [lot . 

Because Eq. [2] contains the terms BW^ {x,m,T , s)^ 
which vary over the Dalitz plot, we can fit the magni- 
tudes aj and the phases 5j of E(^[T]to determine 2/3-1-7 
with only a two- fold ambiguity [71] . Most of the sensitiv- 
ity to 2/3 -|- 7 is expected to come from the interference 
between 6 —> u and 6 — + c transitions leading to D**'^K^ 
final states (with expected r ^ 0.4), and from the inter- 
ference of the former with the b ^ c transition of the 
decay B° ^ D- K*+ . 

The analysis is based on 347 x 10® BB pairs collected 
at the T{4:S) resonance by the BaBAR detector at the 
PEP-II storage rings. The detector is described in de- 
tail elsewhere [lH. In order to estimate signal selection 
efficiencies and to study physics backgrounds, a Monte 
Carlo (MC) simulation based on GEANT4 ^ is used. 

We reconstruct mesons in the decay mode 

K~7T~^7T^. The tracks from £)+ decay are required to 
originate from a common vertex, and the kaon is selected 
using a likelihood based particle identification (PID) al- 
gorithm. The candidates are required to have a 
mass within ±12 MeV/c^ {2a) of the nominal D"*" mass 
[lo| . where a is the experimental resolution. Oppositely 
charged tracks from a common vertex are recognized as 
A'g candidates if they have an invariant mass within ±7 
MeV/c^ (Scr) of the nominal AT^ mass [l3| and a trans- 
verse flight-length significance 4tT greater than zero. The 
TT^ candidate is a track for which the PID is inconsistent 
with its being a kaon or an electron. 

To form B'^ candidates, each £)+ candidate is com- 
bined with a A'g candidate and a 7r~ candidate requiring 
that the three particles originate from a common ver- 
tex. We reject candidates with rn^ (Kgir^) in the win- 
dow [3.40, 3.95] GeVVc^ in order to remove backgrounds 
with non-zero CP content arising from B^ — > D^D^^^ de- 
cays. Using the beam energy in the e"''e~ center-of-mass 
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(CM) frame, two kinematic variables are constructed: the 
beam-energy substituted mass toes = \l 5/4 — P*b^ , and 

the difference between the measured candidate en- 
ergy and the beam energy, = — Here p*g 
and E'^ are the momentum and the energy of the -Brec in 
the CM frame respectively. Candidates with AE in the 
range [—0.1,0.1] GeV and tties in the range [5.24,5.29] 
GeV/c^ are selected. We require [cos 0b |, the absolute 
value of the cosine of the angle between the Bj-oc momen- 
tum and the beam axis, be less than 0.85, and |cos6't|, 
the absolute value of the cosine of the angle between the 
thrust axis of the iJroc decay products and the thrust axis 
of the rest of the event (ROE), be less than 0.95, both in 
the CM frame [H. 

The difference of proper-time At of the two Bs in the 
event is calculated from the measured separation Az be- 
tween the vertices of the i?rcc and the -Stag along the 
beam direction [Q] . We accept events with calculated At 
uncertainty less than 2.5 ps and \At\ < 20 ps. The av- 
erage At resolution is approximately 1.1 ps. The flavor 
of the i?tag is identified from particles that do not be- 
long to the i?rcc using a multivariate algorithm [6|. The 
effective efficiency of the tagging algorithm, defined as 
Q = Ekek{l - 2wkY, is (30.1 ± 0.5)%, where and 
Wk are the efficiency and the mistag probability, respec- 
tively, for each of the six tagging categories k. Untagged 
events contribute to the determination of magnitudes and 
phases of the resonances and they are grouped in a sepa- 
rate seventh category corresponding to the case ^ = in 
Eq. [2]and containing about 38% of the events. 

To further suppress the dominant continuum back- 
ground, which have a more jet-like shape than BB 
events, we use a linear combination T of five variables: 
^0 = Y^iPii -^2 = Si Pi I cos6'j;|'^, the global thrust of the 
event, |cos0t|, and [cos 0b |- Here, pi is the momentum 
and Oi is the angle, with respect to the thrust axis of 
the -Brcc, of the tracks and clusters of the ROE in the 
CM frame. The coefficients of T are chosen to maxi- 
mize the separation between the distributions obtained 
from Monte Carlo simulated signal events and 28 fb^^ 
of continuum events collected at a CM energy 40 MeV 
below that of the T{AS) resonance (off- resonance) , whose 
energy is rescaled to the energy of the beams. The cor- 
relations among the set of measured values of the vari- 
ables (mESi T) arc negligible. Since both T and 
the flavor-tagging utilize the ROE information, the dis- 
tribution of T is correlated with the tagging category. To 
take into account this correlation, we parameterize the T 
distribution for each tagging category separately 

Approximately 7% of selected events contain more 
than one reconstructed signal candidate, arising primar- 
ily from multiple D'^ candidates. We select the one 
having the D-candidate mass nearest to the nominal 
value fio']. For simulated signal events, the entire se- 
lection chain has an efhciency of (9.9 ± 0.1)%, where the 



error is statistical only. 

To separate signal from background and to determine 
their yields, we first perform an unbinned extended maxi- 
mum likelihood fit to the selected on-resonance data sam- 
ple in the variables tties, AE, and !F. The role of this 
first step fit is to extract all the shape parameters, the 
fractions of events by tagging category, and the overall 
yields, which will then be fixed in the subsequent time- 
dependent fit to the Dalitz plot. We define the logarithm 
of the likelihood: 



(4) 



where is the product of the PDFs of mEs, AE, and 
Tk for the event i in the tagging category k. Ntot is 
the total number of events and Njk is the number of 
events of each sample component j: signal (Sig), con- 
tinuum (Cont), combinatoric BB decays (65) and BB 
events that peak in toes but not in the AE signal region 
(Peak). 

The signal is described by a Gaussian function for the 
Toes distribution, two Gaussian functions with common 
mean for the AE distribution, and a Gaussian function 
with different widths on each side of the mean ( "bifur- 
cated Gaussian function") for the T distribution. The 
signal model parameters are obtained from a fit to a 
high-statistics data control sample of B^ — > D^a^ de- 
cays. The selection of these events is similar to signal, 
except that no ifg candidate is required. The decay chain 
— > p^TT^ with — !■ TT^TT^ is reconstructed requiring 
the dipion invariant mass be within ±150 MeV/c^ of the 
nominal p° mass [l3], and the invariant mass of the p 
candidate with the third pion be within ± 250 MeV/c^ 
of the nominal mass [lO|. 

The Toes distributions of the continuum and combina- 
toric BB backgrounds are described by empirical thresh- 
old functions IJ] , while for AE distributions linear func- 
tions are used. The distributions are parameterized by 
a bifurcated Gaussian function for the continuum back- 
ground and a sum of two Gaussian functions for the BB 
combinatoric background. For the latter the parameters 
are determined by BB Monte Carlo simulation. All the 
shape parameters of the continuum background are taken 
from fitting the off-resonance data. The toes distribution 
of the Peak background is parameterized by a Gaussian 
function with the same mean as the signal and a width 
fixed to the value obtained from Monte Carlo simulation. 
The AE distribution is described by an exponential func- 
tion. The T distribution of Peak is described using the 
same PDF as for BB background. 

The yields and the fraction of events for each tag- 
ging category are fitted together with the free shape pa- 
rameters. The yields obtained for each component are 
TVsig = 558 ± 34, A^cont = 13222 ± 226, N^;^ = 5647 ± 213 
and A^poak = 183 ±41, in agreement with the previous 
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result [151. 

The second stage of the analysis is the time-dependent 
fit to the Dalitz plot. For each background component, 
the At distribution is modeled as an exponential, with 
an effective lifetime parameter. To model the detector 
resolution, it is convolved with the sum of three Gaus- 
sians, the sum of two Gaussians and one Gaussian in 
the case of continuum background, BB combinatorial 
background and Peak respectively. The widths of the 
Gaussians, the relative fraction of them, the effective di- 
lution parameters, and the effective lifetimes are deter- 
mined independently from fits to the control samples: the 
off-resonance data sample for the continuum background, 
the BB Monte Carlo sample for BB combinatorial back- 
ground and the Peak component in the BB Monte Carlo 
sample for the Peak. In the case of Peak the lifetime is 
fixed to the B^ lifetime The At parameterizations 
described above for each background component are com- 



signal. The Yskgd parameterization is obtained from a 
weighted average, using the fitted yields, of the shapes 
obtained from the first step fit. 

With the current dataset we are unable to determine 
the magnitudes for the suppressed b ^ u decays. We 



therefore fix the ratio 



0.3 for each res- 



bined in a global time-dependent PDF 77 



±,Bkgd 



obtained 



A{b^c) 

onance in the PDF, which is compatible with the limit 
r < 0.4 (90% C.L.) reported in Ref. fB\. The i?*+(2573) 
magnitude and phase are fixed to the values given in . 
Despite the fact that the b ^ u phases cannot be pre- 
cisely determined they are left free in the fit. All the 
b ^ c magnitudes and phases together with 2/3-1-7 ^'^^ 
free parameters. The whole fitting procedure has been 
validated using high statistic parameterized (toy) Monte 
Carlo samples. 

The fit is performed on events satisfying tties > 5.27 
GeV/c2, |A£;| < 50 MeV and J^- > _2. Results are 
shown in Table [D Figure [1] shows the projections of the 
on-resonance data sample on the two Dalitz plot vari- 



as a weighted average based on thcjitted yields, where ables m^^^o^i) ^^^^ m'^{D^'K^) with the fitted PDFs 



(— ) indicates i?tag — B^ {Bt, 



B°) 



To obtain the PDF describing the Dalitz plot of the 
background in the tagging category k, we use the results 
of the yields fit and calculate for each event a background 
weight [l6j : 



W 1 



S,-Vs.g,jyjfc("^Es,A£;,.Ffc) 



(5) 



where Nj and Yjk are defined as in Eq. [H and Vsig,j is 
the signal row of the covariance matrix of the component 
yields obtained from the likelihood fit. In the absence 
of correlations, M^Ekgd are the background probabilities 
Pskgd/^'totai- Applying these weights to the Dalitz plot 
of on-resonance data we obtain the observed background 
Dalitz plot I?Bkgd- 

For the signal the effect of finite At resolution is de- 
scribed by convolving Eq. [5] with a resolution function 
composed of three Gaussian distributions. Incorrect tag- 
ging dilutes the coefficient of cos(Amc;At) in Eq. [5] by a 
factor (1 — 2'Wi). The parameters of the resolution func- 
tion and those associated with flavor tagging are fixed to 
the values obtained in 6]. 

The expression for the time-dependent Dalitz plot like- 
lihood function is then: 



ln£ = ^ 



ln£+,fc+ ^ ln£_.fe 



BO tag 



S" tag 



(6) 



superimposed. Figure[2K shows the rriES distribution and 
the fitted PDFs for each component, after applying ad- 
ditional requirements on AE and . Besides the value 
of 2/3-1-7, a-i^ important outcome of the analysis is the fit 
of the resonance contributions to the b ^ c part of the 
Dalitz plot. Biases related to the small sample size are 
observed in the measurement of the magnitudes. They 
are estimated using a large number of toy experiments 
generated with the magnitudes values obtained in the fit 
to the on-resonance data sample. 
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FIG. 1: The distributions of m^(A's7r) and m^{D'K) in data 
(solid points). The overall PDF is superimposed. The grey 
full line is the signal PDF, the dashed line is the background 
PDF. 



The likelihood function (>C_.fe) for an event in the 

tagging category k with Stag = B° (i?tac. = -8°) is: 



-'tag — V-^tag 

jk Ti^^ ^n-k ta/c 



C±^k - ^Sig''^±,Sig^Sig + ^Bkgd^Bkgd7'±,Bkgd^Bkgd (7) 

Here Y indicates the product of PDFs for tties, Ai?, and 
^k, 'P±,Sig is the time-dependent Dalitz plot PDF for 



The systematic errors are summarized in Table |TI1 The 
main contribution is related to the parameterization of 
the background Dalitz plot. This effect has been esti- 
mated by repeating the fit with a parameterization ob- 
tained from off-resonance data and BB generic Monte 
Carlo simulation. The systematic uncertainty due to the 
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FIG. 2: a): the tties distribution of on- resonance data (solid 
points) for the global fitted PDF (blue) with the contribution 
of the PDF for each component superimposed: signal (grey 
line), continuum (small dashed line), combinatoric BB (big 
dashed line) and Peak (dotted line). To enhance the signal 
\I^E\ < 0.025 GeV and T > 0.2 have been required, b); 
distribution of the values of 2/3 + 7 fitted on data for different 
hypothesis on the r value, c): variation of the logarithm of 
the likelihood with 2/3 + 7. 



efficiency variation over the Dalitz plot has been evalu- 
ated assuming a flat efficiency. The effect of potential 
CP content of the BB peaking background is taken into 
account assuming the same CP violation structure as in 
the signal with a value = 0.4. The systematic uncer- 
tainties on the signal Dalitz plot come from the variation 
of the r factor (0.3 ± 0.1), of the D*+{2573) magnitude 
(0.02 ± 0.01) and from the introduction of up to 7% of 
a non resonant component. In addition, the masses and 
widths of the resonances have been varied by one stan- 
dard deviation [13] ■ We obtain the systematic uncer- 
tainty arising from imperfect knowledge of the y shape 
parameters and the yields by varying all fixed param- 
eters within their uncertainties. Similar variations are 
applied to the signal and background fractions in each 
tagging category as well as for the At resolution parame- 
ters, the effective lifetimes, the B lifetime and the mixing 
frequency. The systematic uncertainties due to the de- 
pendence of the tagging efficiency on the B flavor, the 
beam spot position and the SVT alignment have been 
obtained following the procedure described in Q . 

Figure[2b shows the dependence of the measurement of 
2/3 -|- 7 on r. For each fixed value of r, a point in the plot 
represents the result of the fit on 2f3+j with its statistical 
error. The error decreases, as expected, for increasing 
r and the central value remains stable. The projection 
on 2/3-1-7 of negative logarithm of the likelihood 
in Figure [2t clearly shows the minimum corresponding 
to the result of the fit and the expected mirror solution 
at +1T rad. Having fixed some magnitudes and strong 
phases, the second solution is disfavored, but it should 



be regarded as equivalent. 

In summary, we present the first results of a time- 
dependent Dalitz plot analysis of the decay i?'^ 
D^K^n^ to determine the Dalitz plot model parame- 
ters and the weak phase 2/3 -|- 7. Assuming r = 0.3 we 
find 2/3 + 7 = (83 ± 53 ± 20)° and (263 ± 53 ± 20)°, where 
the first error is statistical and the second is systematic. 
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Resonance 


Bias correction 
for the magnitude 


Vch magnitude after 
bias correction 


Phase (°) 


K*{892) 

Do (2400) 

7)2(2460) 

7fo*(1430) 
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7f*(1680) 
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-0.011 


1. 
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0. 
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2/3 + 7 


(83 ± 53 ± 20)° and (263 ± 53 ± 20)° 



TABLE IL Sources and sizes of systematic errors. 







A* (2400) 


7)2(2460) 


7^0* (1430) 


A-2*(1430) 


7f*(1680) 


Systematic 


2/3 + 7 




Sc 




Sc 




Sc 




Sc 




Sc 


Bkgd Dalitz plot param. 


16.0° 


0.058 


3.2° 


0.034 


12.1° 


0.088 


9.5° 


0.005 


12.0° 
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10.3° 


Eff. over the Dalitz plot 


5.8° 
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1.9° 
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0.8° 
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19.4° 


CP content of bkgd 
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1.0° 


r 
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2.2° 


0.039 


3.0° 
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0.7° 


0.016 


0.3° 


a(D;+(2573)) 


0.7° 






















m,r 


9.5° 


0.012 


28.0° 


0.011 


6.9° 


0.018 


2.8° 


0.032 


5.9° 


0.036 


9.3° 


y PDF param. 


3.0° 
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1.4° 
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0.4° 


0.007 


0.6° 


0.003 


0.1° 


0.002 


0.5° 


Signal and bkgd frac. 


1.4° 
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2.9° 


0.004 


1.2° 


0.013 


1.4° 


0.008 


0.7° 


0.004 


1.4° 


Yields 


0.1° 


0.003 


1.3° 


0.001 


0.3° 


0.005 


0.4° 


0.002 


0.1° 


0.002 


0.1° 


Tagging and time param. 


2.6° 
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1.4° 


0.001 


0.3° 


0.004 


0.4° 


0.002 


0.2° 


0.002 


0.2° 
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